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The Advantage of Increased Resolution in the Study of Quasar 

Absorption Systems 

Anand Narayanan^, Tor u Misawa^, Jane C. Charlton^, and Rajib Ganguly^ 

ABSTRACT 

We compare a new R = 120, 000 spectrum of PG 1634 + 706 {zqso = 
1.337, my = 14.9) obtained with the HDS instrument on Subaru to a. R = 45, 000 



^ ■ spectrum obtained previously with HIRES/Keck. In the strong Mgll system at 

■ z = 0.9902 and the multiple cloud, weak Mgll system at z = 1.0414, we find 

■ that at the higher resolution, additional components are resolved in a blended 
profile. We find that two single-cloud weak Mgll absorbers were already resolved 

. at -R = 45,000, to have b = 2-4 km s~^. The narrowest line that we measure 

in the R = 120,000 spectrum is a component of the Galactic Nal absorption, 
(yQ ■ with b = 0.90 ±0.20 km s~^. We discuss expectations of similarly narrow lines in 

^ . various apphcations, including studies of DLAs, the Mgl phases of strong Mgll 

O ■ absorbers, and high velocity clouds. By applying Voigt profile fitting to synthetic 

lines, we compare the consistency with which line profile parameters can be accu- 
rately recovered at i? = 45, 000 and R = 120, 000. We estimate the improvement 
gained from superhigh resolution in resolving narrowly separated velocity com- 
. ponents in absorption profiles. We also explore the infiuence of isotope line shifts 

and hyperfine splitting in measurements of line profile parameters, and the spec- 



tral resolution needed to identify these effects. Super high resolution spectra of 
quasars, which will be routinely possible with 20-meter class telescopes, will lead 
to greater sensitivity for absorption line surveys, and to determination of more 
accurate physical conditions for cold phases of gas in various environments. 



1. INTRODUCTION 

With the advent of high resolution spectroscopy on 8-meter class telescopes in the early 
1990's came a revelation in the study of quasar absorption line systems. Previous to this time, 
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several large surveys were conducted with a typical resolution of R — A/ A A = 1000-3000, 
corresponding to 100-300 km s~^. These surveys tabulated the redshift path densities of 
Mgll, Civ, and Lja lines over the redshift range accessible from optical telescopes (Lanzetta 
et al. 1987; Sargent et al. 1988; Caulet 1989). However, most individual systems were not 
resolved, so that the kinematics of the different metal lines could not be compared system 
by system. It is not surprising that an increase of resolution to = 45, 000 (~ 7 km s^^) 
provided detailed views of the multiphase structure of gas inside and around galaxies, since 
parcels of gas within the different components of galaxies are often moving at several to 
tens of km s~^ relative to one another. The increase of resolution from thousands to tens 
of thousands also brought a dramatic increase in the sensitivity of surveys. This allowed 
measurements of the mctallicity in the most diffuse gas sampled by weak Lya forest lines 
(e.g., Cowie et al. 1995), and for deeper surveys of metal lines as well (e.g., Churchill et al. 
1999; Misawa et al. 2002). 

Presently, hundreds of quasar spectra have been obtained with R ~ 45, 000 using the 
High Resolution Echelle Spectrograph on Keck I (HIRES, Vogt et al. 1994), the High 
Dispersion Spectrograph on Subaru (HDS, Noguchi et al. 1999) and the Ultraviolet-Visual 
Echelle Spectrograph on the VLT (UVES, Dekker et al. 2000). There is undoubtedly 
much to be gained from detailed studies of the absorption hnes detected in these databases. 
However, in this paper we look toward the future and ask what advances we should expect in 
studies of quasar absorption lines with an increase in resolution from the presently attainable 
value of 7 km s~^ to 0.5-2 km s~^, which is within the reach of 20-meter class telescopes 
for hundreds of quasars. This follows a study by Tappe & Black (2004) of the quasar 
PKS 2145 + 067 at i? ~ 100, 000 with UVES on the VLT UT2 telescope. That study focused 
on the z = 0.79089 strong Mgll absorber, and measured b parameters in the range 1- 
7 km s~^, showing considerably more structure than previous spectra (Churchill et al. 2001) 
with R ~ 45, 000. More recently, Chand et al. (2006) considered the effect of observing at 
R > 100, 000 for studies of variation in the fine-structure constant. 

In addition to the increased sensitivity to weak lines, gained at higher resolution, there is 
reason to expect we should see substructure on km scales. Distinct small-scale structures 
are seen in the interstellar medium of the Milky Way on these velocity scales (Points et al. 
2004). Very narrow components, indicating cold clouds, are seen through 21-cm absorption in 
some damped Lya systems (Lane et al. 2000). Also, detailed models of transitions detected 
in DLA systems indicate the presence of a cold phase (Wolfe et al. 2003). The large Mg 
to Mg+ ratios in many components of strong Mgll absorbers, most of which are produced 
by lines of sight through L* galaxies, could be explained by a cold phase that produces the 
bulk of the Mgl absorption (Churchill et al. 2003). In the case of the z = 0.9902 system 
toward PG1634 -|- 706, Ding et al. (2003) proposed cold clouds with Doppler parameters of 
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b{Mg) ~ 0.75 km to explain the strong Mgl absorption, but these could not be resolved 
at FWHM~ 7 km s~^. Finally, VLA observations of Milky Way high velocity clouds indicate 
that they have T < 900 K, and most Milky Way lines of sight show narrow Call lines from 
HVCs (Richter et al. 2005). However, at FWHM~ 7 km s~^ these Call lines were unresolved. 
Richter et al. (2005) suggested that if the Milky Way is typical of spiral galaxy halos these 
narrow absorption features should be commonly seen in the profiles of strong Mgll absorbers 
and Lyman limit systems. 

For the brightest quasars it is already feasible to use the highest resolution settings on 
8- meter class telescopes (~ 40 with my < 16 for non-exorbitant exposure times) to study 
quasar absorption line systems. In this paper, we present an exploration of the power of 
it! = 120, 000 in studying four Mgll absorbers along the line of sight toward the magnitude 
mv — 14.9 quasar PG1634 + 706. Our observations were carried out in the service mode 
with the HDS on the Subaru Telescope. 

The paper begins, in § 2, with a general consideration of how accurately line parameters 
can be determined at super-high resolution. In § 3, we present our Subaru/HDS observations 
of several absorption line systems in the spectrum of one of the brightest quasars, PG 1634-1- 
706, and discuss the implications of our R = 120, 000 observations of these systems. We also 
show the Milky Way absorption along this line of sight. In § 4, we speculate about the 
progress that will be made in the study of quasar absorption lines when many quasars are 
observed at R> 100, 000. 

2. SIMULATION OF VERY NARROW ABSORPTION LINES 

In gas with temperatures T < 8000 K, most metals would produce absorption features 
that are narrower than btherm — 2 km s""^. Figure 1 illustrates progressively narrow absorp- 
tion profiles of MgllA 2796 hne with column density A'"(Mg+) = 10^^ cm~^ and with Doppler 
parameters &(Mg+) from 9 km s~^ to 0.01 km s~^ at high {R — 45, 000) and super-high spec- 
tral resolutions (R = 120,000). The various profiles were produced by convolving a model 
absorption feature with Gaussian kernels of FWHM = 6.6 km s^^ and FWHM = 2.5 km s~^, 
to simulate the resolving powers corresponding to the two spectral resolutions. Comparing 
the profile shapes, it is apparent that, even at high resolution, the width of intrinsically 
very narrow lines is often misrepresented. A measurement of the b would result in a value 
significantly higher than the true value, which further affects the estimation of the upper 
limit for the gas temperature. To explore the effect of resolution further, we evaluate the 
consistency with which Voigt profile parameters (i.e., column density and Doppler parame- 
ter) can be recovered from very narrow absorption profiles at two resolutions; R — 120, 000 
and it! = 45, 000 in synthetically simulated spectra. Both these parameters are important for 
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deriving ionization structure, gas phase metallicities, and other physical conditions for the 
absorber. In addition, using this approach, we also consider two other effects: (1) blending 
of lines, and (2) isotope shifts and hyperfine splitting in absorption lines. We evaluate the 
advantage of increased resolution in recognizing these effects and also in determining how 
they influence Doppler parameter measurements of absorption profiles. 



For our simulations wc synthesized Mgll A 2796 lines with two sets of values: (1) column 
density of A^(Mg"'") = 10^^ cm^^ and 6(Mg+) = 2.0 km s~^, corresponding to a rest-frame 
equivalent width hmit of W^(2796) = 0.03 A, and (2) N{Mg+) = 10^^ ^.^^-2 6(Mg+) 
= 1.0 km s~^, corresponding to 1^^(2796) = 0.02 A. The former characterizes the typically 
measured profile parameter for the unique class of weak metal line absorbers known as 
"single-cloud weak Mgll" systems (Churchill et al. 1999; Rigby et al. 2002; Charlton et al. 
2003). The latter would be representative of even colder phases of gas (T < 1460 K) at low 
ionization. Our choice of these values for 6(Mg+) for the simulations was also motivated 
by a desire to find the minimum Doppler width at which we expected to resolve lines at 
R = 45, 000, and to show the improvement with R = 120, 000. Thus the simulations are 
driven by the characteristics of the particular instruments that we compare. 

In generating a synthetic spectrum, an ideal absorption line for the chosen parameters 
was first simulated using a Voigt function to represent the dependence of absorption on 
wavelength. The synthetic profile was initially oversampled to assure accurate representation 
of narrow lines. It was then convolved with a Gaussian instrument kernel, as described above. 
The convolved profile was then resampled at a rate of three pixels per resolution element to 
enable comparison with observed spectra for a particular R. Finally, Poisson noise was added 
to the spectra, corresponding to a specified S/N per pixel. To determine the parameters N 
and b we performed Voigt Profile fits on the absorption lines in 10, 000 simulated spectra. 
Figure 2 shows a sample from those simulated lines. The fit parameters were derived by 
first generating an initial model profile using an automated fitting routine {AUTOVP, Dave 
et al. 1996), and then adjusting parameters to minimize the chi-square of that initial model 
{MINFIT, Churchill 1997), after convolving with the instrumental spread function. 

Since we want to evaluate the balance between S/N and resolution, we consider pairs of 
equal length exposures with R — 120, 000 and R — 45, 000. The S/N ratios were calculated 
using the relation: 



2.1. Very Narrow Single Component Lines 
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where F is the expected number of photons incident on the spectrograph sht, 5 A is the 
grating dispersion (A/pixel), and f{s) is the sht throughput for the assumed seeing (s). For 
simphcity we assume that the incident numbers of photons are the same (i.e, F45 = -^120)5 
though for real systems this is a function of the telescope throughput and collecting area. (We 
note also that this scaling assumes object-limited observations with negligible contributions 
from the background and detector read-noise.) We assume that the grating dispersions 
are 0.04A and 0.016A for the R=45,000, and 120,000 spectra, respectively. To calculate 
the sht throughput, we consider a seeing of 0.6" (FWHM, with a Gaussian PSF) with 
slit widths of 0.861" and 0.3" for the R=45,000, and 120,000 cases. These slit widths 
were chosen because they represent the two instrumental set-ups used in our observations of 
PG 1634 + 706 with Keck I/HIRES and Subaru/HDS. In our R = 120,000 observation of 
PG 1634 + 706 using Subaru/HDS we were able to obtain a S/N ~ 55 pixel"^ for a 1 hour 
exposure under 0.6" seeing condition. Using the relation described above, this scales to a 
S/N of 128 pixel^^ for a spectral resolution oi R = 45, 000 for identical exposure time and 
seeing condition. 

Figure 3 presents the distributions of the measured values of b and for 10, 000 re- 
alizations of simulated spectra with N = 10^^ cm~^ and b = 2.0 km s~^ lines for various 
combinations of R and S/N. Doppler parameter characterizes the width of the Voigt profile 
which in turn influences the gas temperature estimation. Therefore we compare the consis- 
tency in the measured 6- values for each of the 10, 000 simulated spectra. For the it! = 45, 000 
simulations, we flnd the measured fo-values are distributed with a large tail to higher values. 
The mode of this distribution also falls at a value that is higher than the true 6-parameter 
and it shifts to even higher values as the S/N ratio is decreased. In a small subset of cases 
the Voigt profile models from our fitting procedure yielded results that have large formal 
uncertainties, a (6), in the measurement of b. We find that measurements with a{b) > b 
are poorly constrained, often with the data in one or more pixels deviating by more than 
3 a from the model fit. These models are shown in the distribution plots of Figure 3 as a 
separate category. For the S/N — 80 case, the small peak at ~ 0.95 km s~^ results from 
a pixelization effect. Typical line profiles with the appropriate sampling are detected in a 
small discrete number of pixels (~ 1-2), which can lead to bias toward particular values of 
b. For these profiles the noise also sometimes conspires to produce a b value even narrower 
than the instrumental spread function. In such a case, the true b value is likely to be small, 
but the measurement cannot be claimed as accurate. These simulations therefore reveal the 
inadequacy of R = 45, 000 data in allowing the recovery of a narrow 6-value. 

In comparision, for the super-high resolution synthetic spectra, the distribution of the 
measured b is significantly narrower and the mode is very close to the true value. Additionally, 
the intrinsic value remains the most frequently recovered 6-parameter even as the S/N ratio 
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is reduced. For the S/N that are compared, in almost all 10, 000 cases our fitting procedure 
gave a well constrained measurement of the 6-value. Table 1 lists, for the various simulations 
with the two resolutions, the median value of b, with asymmetric 1-a error bars, and the 
fraction of Voigt Profile fit models with accurate fits [b > crib)]. We conclude that at 
R = 120, 000, a line profile with a true 6-value of 2 km s~^ is very accurately recovered. At 
i? = 45, 000, there is a significantly wider distribution of measurements, with a non-negligible 
chance of obtaining b significantly greater than 2 km s~^, and even occasional measurement 
of 6 ~ 1 km s"^ due to pixehzation effects. 

Figure 4 shows the distribution of the column density measurements for the same set 
of simulations. The distributions are broader for the R — 45, 000 case, with mode values 
smaller than the true value. Measurements thus preserve equivalent width in the sense that 
a line with an N measurement smaller than the true value will have a measured b larger than 
the true value. The Voigt profile models that gave large a{b) also yield large uncertainties 
in column density a (N). Table 1 also summarizes our measurements of iV- values for the 
various simulations. Comparing values in Table 1, it is evident that the constraints from 
these model fits arc poor. Wc conclude that there is a benefit of observing at i? = 120, 000 
for lines with 6 = 2 km s^^, however it is likely that even at i? = 45, 000, the derived value 
will be roughly correct. 

The real advantage of higher resolution is for even narrower lines. We performed a 
similar series of simulations for Mgll A2796 lines with N — 10^^ cm~^ and b — 1 km s~^. 
The results for the distributions of b parameters for various S/N values a,t R — 120, 000 
and R — 4:5, 000 are presented in Figure 5. The large spread in the measured value of b for 
= 45, 000 is indicative that this spectral resolution is inadequate to faithfully recover fine 
widths of narrow line profiles. We again, for R = 45, 000 find a large number of measurements 
with large o"(6) values, for measured 6 of ~ 0.95 km s~^. This is again explained by a 
pixehzation effect, but in this case the measurement happens to be very close to the true 
value. Even if we had chosen a true value of b between 1 and 2 km s^^, this peak would still 
be at 6 ~ 0.95 km s^^. Nonetheless, a measurement of a small b value does imply a true 
narrow line-width, even if the measured value is biased by the pixehzation. 

Even for the R = 120, 000 spectra, the mode of the distribution for measured-5 is 
at a value slightly larger than the true b although in general the measurements are well 
constrained. Table 1 lists the median values of Doppler parameters and column densities, 
and their 1-a errors, for various S/N values at the two resolutions. The distributions of 
A'^-values are compared in Figure 6. 

The significance of increased resolution is clearly seen from Figures 3 through 6. The 
determination of b is also sensitive to the sampling, where b greater than roughly a pixel 
is required. In comparison, the precision in measurements from lower resolutions can be 
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improved only by using substantially higher S/N spectra thus requiring very long exposure 
time. Furthermore, we find that similar equivalent width limits are attained in the R = 
120, 000 and R = 45, 000 spectrum for equal exposure lengths though the S/N ratio of the 
former is substantially lower. For example, the 5o" equivalent width limit for R = 45, 000, 
S/N = 80 and R = 120, 000, S/N = 35 is 5 mA/pixel. Also, as we discuss in the following; 
sub-section and in § 3, observations with R = 120, 000, corresponding to a resolution of 
2.5 km s~^, are very important for resolving lines and separating velocity sub-structures. 

2.2. Blending of Lines 

To probe the nature of small scale structure in absorbers it is essential to resolve closely 
separated velocity components. Here we focus on making estimates of the scale at which 
superhigh resolution offers a real improvement over intermediate resolution spectra. We 
chose a simple scenario of two absorbing clouds, with different physical properties, marginally 
separated in velocity. The column densities and Doppler parameters of the simulated lines 
were selected to match closely with those of the weak multiple cloud system at z = 1.0414 
in the spectrum of PG 1634 + 706, so that our choice is not completely arbitrary. This 
multiple cloud system is discussed in § 3.3. The two components have column densities 
A^(Mg+) of 10^^ and 10^^'^ cm^^ and Doppler parameters 6(Mg+) of 4.0 and 3.0 km s^^, 
respectively. This corresponds to an equivalent width ratio of approximately 3 : 1 between 
the two components. In a set of 1000 simulations with a, S/N oi 80 pixel"^ a,t R — 45, 000, 
and with the two clouds separated in velocity by Av=10 km s~^, we found that in ~ 84% of 
the cases Voigt Profile fitting was able to distinguish the two components in the absorption 
profile. For the same S/N ratio, with the velocity separation reduced to Av=7km s~^, 
in ~ 95% of the cases the simulated line was fit with just a single broad Voigt profile 
with a mean value of (6(Mg+))= 6.3 km s^^. This translates to a larger upper limit for 
the estimated temperature of the gas (T < 58000 K) compared to the true value of T 
< 20000 K. For fines with the same intrinsic properties, but simulated at i? = 120, 000 
and S/N — 35 pixel"^ (corresponding to an identical exposure time), we found that Voigt 
profile fits recovered both kinematic components in ~ 97% of the cases with a velocity 
separation of Av=7 km s"^, and in 100% of the cases with the velocity separation increased 
to Av=10 km s~^. Furthermore, increasing the S/N ratio of the intermediate resolution 
spectra did not compensate for the lower resolution. For example, when the two components 
were separated in velocity by 7 km s~^, even with S/N = 160 pixel"^ at R = 45,000, the 
kinematically separate components were resolved only in ~ 40% of the cases. Figure 7 shows 
a sample of the simulated lines at the two resolutions. From these results we infer that even 
at fairly high S/N ratios, intermediate resolutions can be completely inadequate in revealing 
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closely blended velocity structures on scales smaller than 8 km s~^. Further, such unresolved 
lines can lead to overestimation of the temperature of the absorbing gas. More generally, we 
conclude that, if two lines are narrower than their velocity separation, they can be resolved 
by a Voigt profile fit if the resolution is less than ~ 70% of the velocity separation. The exact 
criterion for separating blends is dependant on the quality of the spectra and the Doppler 
parameters of the two lines. 

2.3. Contribution from Isotopes and Hyperfine Splitting 

Magnesium has three isotopes with mass numbers 24, 25 and 26 a.m.u, with relative 
abundances of 78.99%, 10.00%, 11.01%, respectively. The wavelengths of the resonance fines 
for these isotopes differ from one another. Furthermore, ^^Mg shows hyperfine splitting. The 
relative shifts in rest frame wavelength for each line of the resonance doublet are measured 
to be A24 - A+ = 0.0014 A, A24 - X25 = 0-0061 A, and A24 - A26 = 0.008 A (Morton 2003). 
In velocity space this corresponds to 0.1500 km s~^, 0.6540 km s~^ and 0.8577 km s~^, 
respectively. For such narrow velocity separations, the effect of isotope shifts and hyperfine 
splitting becomes distinguishable only at extremely high spectral rcsohitions {R > 500, 000), 
as illustrated in Figure 8. To estimate the scale at which this effect is likely to inffucncc our b- 
value measurements, we simulated very narrow synthetic lines with true 6(^^Mg+) of 0.3, 0.7, 
and 1.5 km s~^. The 6-values of the other isotope fines were thermally scaled. The simulated 
fines had S/N — 35 pixel"^ at R — 120, 000. Figure 9 shows the dispersion in the incidence of 
measured b for the various cases. Of the 10, 000 simulated fines with 6(^^Mg+) = 0.3 km s""^, 
only ~ 10% of the measurements yielded a 6- value less than 0.5 km s~^, with the recovered 
6-value almost always larger than the true value. For a higher true &-value of 0.7 km s~^, the 
measurements were still distributed across a range from 0.7 to 1.5 km s~^. The peak of this 
distribution, however, was closer to the true-value as compared to the previous case. For 
the final case, with b = 1.5 km s~\ the spectral resolution became adequate to repeatedly 
measure values close to the true b value along with consistently yielding measurements that 
are refiable (cr(6) < b). In this case, although the isotope and hyperfine fines were still 
contributing to the final absorption profile, their contribution did not induce a significant 
change in the line width as the individual line widths were now larger than the largest 
separation between the isotope lines. 

With these results we infer that isotopic fine shifts and hyperfine spfitting are to be 
considered as influencing our measurements only when the true line width is really narrow 
(T < 700 K, 6(Mg"'") < 0.7 km s~^). Furthermore, for such narrow Mgll lines, identifying 
the effect would require spectral resolutions oi R> 500, 000. 
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3. SUPER HIGH RESOLUTION OBSERVATION OF PG 1634 + 706 

Wc obtained a super high resolution spectrum of PG 1634+706 with the High Dispersion 
Spectrograph (Noguchi et al. 1999, HDS) on Subaru. Using a 0.3" slit width, we were able 
to attain R— 120, 000 spectrum with a sampling rate of 3 pixels per resolution element. Our 
1 hr observation yielded a spectrum with an average S/N 30 pixel"^ over the entire spectral 
range. The reduced ID spectrum was extracted in a standard manner with the IRAF software 
^. Wavelength calibration was done using a standard Th-Ar spectrum. The spectrum 
was continuum fitted with a third-order cubic spline function and normalized. Previously 
identified weak and strong Mgll absorption systems in a i? = 45, 000 Keck/HIRES spectrum 
of PG 1634 + 706 made this observation an ideal test case. The HIRES data aquisition 
and reduction procedure is described in (Churchill et al. 2001). The short exposure on 
HDS reached an equivalent width hmit several times smaller than our exposure of twice its 
length with the HIRES instrument at a lower resolution of = 45, 000. Details of the two 
observations are listed in Table 2. Figures lOa-lOc show the wavelength regions that are 
covered in both the HDS and HIRES spectrum. 

In the following subsections, we compare detections of the strong, single cloud weak , 
and multiple cloud weak Mgll absorption features in the super-high resolution HDS and the 
high resolution HIRES spectra for this quasar. We also consider the Milky Way absorption 
along this line of sight. We discuss the particular case of these systems, but also generalize 
to discuss how studies of the same type of system would be affected by higher resolution 
observations. Since the scope of this paper is to compare observations taken at different spec- 
tral resolutions, we defer detailed analyses of the absorption-line systems presented herein 
to future work. 

3.1. Strong Mgll absorber at z = 0.9902 

The z — 0.9902 system along the line of sight to the quasar PG 1634-1-706 has previously 
been studied using a combination of the R — 45, 000 Keck/HIRES spectrum, data from 
HST/FOS (Churchill et al. 2000; Charlton et al. 2000) and from HST/STIS (Ding et al. 
2003). Churchill et al. (2003) found that this strong Mgll doublet could be adequately fit 
with five blended components. With the higher resolving power of HDS, additional kinematic 
components were revealed in the Mgll absorption feature. Our Voigt profile fit decomposes 
the Mgll AA2796, 2803 profile into seven components. Figure 11 shows the HIRES and HDS 
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data, comparing differences in the profile shapes for the two resolutions. The asymmetrical 
shapes for the MgllA2796 line at Av ~ 22 km and Av ~ 8 km from the flux-weighted 
system center in HIRES spectrum were indicative of narrowly blended features which are 
clearly distinguished in the profiles in the HDS spectrum. Table 3 presents our Voigt fitting 
parameters for Mgll AA2796, 2803 from the Keck/HIRES and the Subaru/HDS observations. 
This ability of superhigh resolution to resolve additional components in strong metal line 
absorbers was also noted by Chand et al. (2006). 

Strong Mgll absorption systems usually have absorption from gas in different phases 
(different densities and temperatures) residing within an impact parameter of 35h~^ kpc of a 
luminous galaxy with L > 0.05 L* (Bergeron & Boisse 1991; Steidel et al. 1997). Dissecting 
the Mgll absorption profile to reveal kinematic structures at small velocity separations of 
Av < 1 km is therefore useful in developing a physical understanding of the spatial 
distribution of gas in galaxies along the line of sight. If the observed Mgll absorption 
signatures are characteristic of disk ISM (Lanzetta & Bowen 1992; Charlton & Churchill 
1998), then measuring structure on the smallest scales is important to determine if the 
interstellar medium of high rcdshift galaxies follows a complex fractal structure similar to 
Milky Way (e.g., Elmegreen 1997). Resolving closely blended components, as in the z ~ 
0.9902 system is also essential in developing ionization models that provide strong constraints 
on physical parameters. This is especially relevant when there is a difference or gradient in 
metallicity, abundance pattern or ionization from cloud to cloud. 

Churchill et al. (2003) studied 23 Mgll systems at z ~ 1 and found that many strong 
Mgll systems have observed column density ratios of N(Mg)/N(Mg+) > 0.01. This value 
is the maximum expected if the absorption is taking place in gas at a single phase under 
photoionization equilibrium. In developing a model for this z = 0.9902 absorption system 
based on the R = 45, 000 Kcck/HIRES observations. Ding et al. (2003) found that the 
Mgl absorption line in this system is too strong to arise in the same phase as Mgll for all 
possible ionization parameters. To reconcile this discrepancy. Ding et al. (2003) invoked 
an additional lower ionization phase that gave rise to the majority of Mgl absorption, but 
which produced small amounts of Mgll. The lower ionization phase also produced most of 
the neutral hydrogen absorption from this system, with a range 17.3 < logA^(H) < 17.7 
for the five Mgl clouds. In this model, the Mgl phase gas has a very low temperature 
(T ~ 800 K, 6 ~ 0.75 km s'^), high density (uh ~ 200 cm-^) and small size (~ 100 AU), 
reminiscent of a cold pocket embedded in a warm, ionized inter-cloud medium as in the 
ISM of galaxies. The spectral resolution of HIRES was inadequate to show the narrow line 
widths that we might expect from such cold, small-scale clumps of gas. However, the present 
Subaru/HDS spectrum provides a direct test of the Ding et al. (2003) prediction for this 
z = 0.9902 absorber. 
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Figure 12 shows the absorption from Mgl A2853 in the z — 0.9902 system, as seen in 
the HDS spectrum. The Voigt profile fit constraints on the various hne features is hsted in 
Table 3. For the three clouds used to fit the Mgl profile, the line widths are not as small as 
what Ding et al. (2003) had predicted in their models (see their Table 1). For the line width 
predicted by their model, the Mgl absorbption has to be from gas at temperature T < 800 K. 
Our simulation results from § 2.3 show that in instances where the true-6 value is very small, 
Voigt profile measurements can yield higher values because of unresolved isotopic hues. To 
investigate if our measurements of 6- value for the various components of Mgl is influenced by 
this effect, we resort to simulations. The isotopic lines of Mgl occur at A24 = 2852.9635 A, 
A25 = 2852.9655 A and A26 = 2852.9674 A (Beverini et al. 1990). The Mgl 2853 A hne does 
not have hyperfinc structure. In our simulations, we used the profile parameters predicted by 
the model (N(Mg) = 10"'^^'"'^ cm^^, 6(Mg) = 0.75 km s^"*^) for the isotopic line of ^^Mg, with the 
6-parameters of the other two isotopic lines thermally scaled. The spectra were simulated for 
a S'/A^ = 35 pixel^^ at i? = 120, 000, with the lines redshifted to match the redshift of the 
system. Figure 12 shows the distribution in the incidence of measured h from a simulation 
of 1000 hues. It can easily be noticed that, although the measured h is almost always larger 
than the true h of 0.75 km s~^, a measurement of 2 km s~^ or higher (as in the observed 
HDS spectrum) never occurs. It is therefore to be concluded that the hypothesis of a two 
phase structure as an explanation of the Mg to Mg+ ratio is therefore not consistent with 
the high resolution observations. It is clear from § 2 that a component with h = 0.75 km s~^ 
will not be measured to have b as high as several km s~^ at i? = 120, 000. Therefore, another 
explanation must be sought for this system. Tappe & Black (2004) suggested that a UV flux 
contribution from the absorbing galaxy and/or Mg + charge transfer could significantly 
influence the Mg to Mg+ ratio. Their R = 100, 000 VLT/UVES observation of a different 
system, at z = 0.79089 toward PKS 2145 + 067, yield 6(Mg+) = 1.02 ± 1.65 km s"^ for the 
one component detected in Mgl in a spectral region with S/N 30. The error is too large 
to draw a conclusion in that case. Another possibihty for the z — 0.9902 system is that there 
are more components in the Mgl than we are able to resolve. 

We conclude that our R — 120, 000 spectrum does not provide support for a cold phase 
model for the z = 0.9902 system, though some modified version could apply. Direct tests 
with super high resolution for other similar systems are needed. One excellent candidate 
is the z = 0.4523 multiple-cloud, weak Mgll system toward the quasar HE 0001 — 2340 
{z = 2.28, my = 16.7), which has detected Fcl, and Cal as well as Mgl (Jones et al. 
2006). Cold phase models may also have significant consequences for the study of damped 
Lya absorbers (DLAs) and sub-DLAs. Lane et al. (2000) have resolved narrow components 
{b ~ 2-3 km s^^) in an Hi 21cm study of the DLA system at z = 0.0912 toward B 0738+313. 
Mgl lines in these systems are expected, from thermal scaling, to be significantly narrower. 
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Also, more recently, Srianand et al. (2005) report measuring kinetic temperature of ~ 150 K 
from H2 lines in three DLAs. Such a low gas phase temperature can produce narrow metal 
lines, with Doppler width of ~ 3 km s^^, that can be accurately recovered with superhigh 
resolution spectra. A fair fraction of DLAs have detected molecules (e.g., Petitjean et al. 
2000) , indicating directly that they have a cold phase. Other DLAs and Lyman limit systems 
may also be produced by lines of sight through cold regions surrounding molecular clouds, 
or regions that did not quite reach the threshold for star formation. 

3.2. Single-Cloud Weak Mgii absorbers 

Most single-cloud, weak Mgll absorbers (rest-frame equivalent width Wr(2796) < 0.3 A) 
are not associated with luminous galaxies, although models designed to fit the Lyman series 
line and the observed lack of Lyman break constrain the metallicity of the absorbing gas to 
be at least 10% solar. Photoionization models suggest that the structures responsible for 
weak Mgll absorption might be unstable over astronomical timescales, because of pressure 
imbalance between a high density, low ionization {p/k = 10^^ cm^"^ K^^) and a low density 
high ionization {p/k = 10^^ cm^^ K^^) phase (Charlton et al. 2003; Narayanan et al. 2005). 
The conclusions from photoionization models arc based on the observed values of column 
density and Doppler parameter for the low and high ionization transitions. The best con- 
straints on these parameters were mostly from high resolution data at ~ 45, 000 or lower. 
The b parameters measured at ~ 45, 000 (6 ~ 2 or 3 km s~^) allowed for temperatures 
of thousands of Kelvin, similar to the temperatures inferred for the high ionization phases. 
If the true Doppler parameter is smaller than that derived from the HIRES measurements, 
then it would imply that the low ionization phase has a much lower temperature than the 
high ionization phase. This influences the inference about their stability. Line widths nar- 
rower than ~ 2 km s~^ (T < 15, 000 K) cannot be accurately measured at i? = 45, 000. 
Therefore the temperature estimations for the low ionization phase in these absorbers may 
not be accurate. Our R — 120, 000 Subaru/HDS spectrum allows us to test the assumption 
that these single-cloud weak Mgll hues are resolved. 

For the two single cloud weak Mgii absorbers ai z = 0.9056 and z = 0.8181 found in 
the PG 1634 -|- 706 spectrum, we find from the HDS spectrum that the weak Mgll hnes are 
not any narrower than values measured at R 45, 000 (Figures 13). The column densities and 
Doppler parameters from simultaneous Voigt profile fits to Mgll AA2796, 2803 are listed in 
Table 4. Thus at least for these cases, the lines were already sufficiently resolved at R = 
45, 000 to allow measurements of the 6- values that were accurate within the errors. These 
measurements are in concordance with the instability inference based on photoionization 
model. Analogs of weak Mgll absorbers are known to exist in the present day universe 
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( Narayanan et al. 2005) as well as at intermediate ~ 1, Churchill et al. 1999) and 
higher redshifts {z ~ 2, Lynch et al. 2006). Therefore confirming the transient nature 
of these systems is significant in addressing the question of the evolution of the processes 
that created these absorbing structures. A larger sample of weak Mgll absorbers should be 
measured with super high resolution in order to verify that we can generalize this conclusion. 
However, the real value in R > 100, 000 observations in the study of weak Mgll absorbers 
will be the enhanced sensitivity which will allow us to estabhsh the shape of the distribution 
function for the weakest absorbers, and to infer their physical properties. 

3.3. Multiple-Cloud Weak Mgii absorbers 

Approximately 35% of weak Mgll absorbers have multiple clouds (Rigby et al. 2002). 
Based on their kinematics, these multiple-cloud, weak Mgll absorbers can be described as 
"kinematically compact" or "kinematically spread". The kinematically compact systems 
have been hypothesized to arise in dwarf galaxies, and the kinematically spread systems in 
the outskirts of spiral galaxies (Masiero et al. 2005; Ding et al. 2005), but it is likely that the 
situation is more complicated. Varying degrees of line blending in the absorption features 
implies that the Voigt profile fit results arc sometimes going to be sensitive to resolution. 
This is illustrated by the z = 1.0414 multiple-cloud, weak Mgll system in the spectrum of 
PG 1634-1-706. Based on HIRES data, this weak Mgll absorber was classified as a multiple 
cloud weak Mgll system with 4 discrete Mgll clouds (Churchill et al. 2003). Two additional 
weak components were resolved in Voigt profile fits to the higher resolution data in the same 
region, as shown in Figure 14. Zonak et al. (2004) produced models for the system based 
on the HIRES data along with R = 30,000 ultraviolet spectra obtained with HST/STIS. 
It is apparent from our HDS spectrum, in Figure 15, that there are two subsystems in 
Mgll, separated by ~ 150 km s^^. For the blueward subsystem (see Zonak et al. 2004), 
Mgll A2796 was not covered in the HIRES spectrum, and Mgll A2803 was not detected, but 
several other metal-fine transitions (e.g. Sim 1207, Si IV AA1394, 1403 etc) were detected in 
the HST/STIS spectrum. The greater sensitivity of the R = 120, 000 HDS spectrum led to 
detection of Mgll A2803 in the blueward subsystem, which was a prediction of the models 
of Zonak et al. (2004) 

Enhanced resolution allows the detailed shapes of absorption profiles to be surveyed, 
and blends to be separated. Some single-cloud weak Mgll absorbers appear asymmetric 
when observed at it! = 45, 000 (Churchill et al. 1999), suggesting they may have more than 
one Voigt profile component. More generally, it is of interest what fraction of the single- 
cloud weak Mgll absorbers are truly fit by single Voigt profile components when viewed 
at higher resolution. The profiles of the z = 0.9056 and z = 0.8181 systems observed at 
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R — 120, 000 here, still are fit adequately by a single component, suggesting a very simple 
structure. Observing some asymmetric single-cloud weak Mgll may reveal some to be more 
similar to kinematically compact multiple cloud weak Mgll absorbers, or they could just be 
variations of ordinary single-cloud weak Mgll absorbers. Determinations of their metallicities 
and comparisons of kinematics of other transitions should allow a separation of these classes. 

3.4. Milky Way Absorption and High Velocity Clouds 

In the HDS spectrum, Galactic Nal is detected, as shown in Figure 16. Voigt profile 
fitting parameters are listed in Table 5. For the narrowest component, at vlsr = —32 km s~^ 
^ we measure 6 = 0.9 ± 0.2 km s~^. This shows that such narrow fines as 6 ~ 1 km s~^ can 
be measured a,t R — 120, 000. At the same time, we also note that the hyperfine splitting in 
Nal lines are likely to affect our measurement of such a narrow line width. The difference 
in rest-frame wavelength from the basic hyperfine splitting of the Dl and D2 lines has been 
measured to be 0.0213 A and 0.0198 A, respectively (Morton 2003). This corresponds to 
shifts in velocity of 1.0827 km s^^ and 1.0075 km s^^, respectively. As these values are 
very close to our measurement of 6 ~ 1 km s^^, we expect the truc-6 to be lower than this 
measured vahic. However, observing this effect by resolving the line would require spectral 
resolutions that arc significantly larger than our HDS measurement at i? = 120, 000 (as 
discussed in § 2.2.) 

The PG 1634 + 706 line of sight passes through the well-studied high velocity cloud 
complex C. Upon inspection of a HST/STIS E230M spectrum, we find absorption detected 
in Mgll, Mgl, and Fell at vlsr ~ -115 and -155 km s'^ We do not detect NalAA5892,5898 
at these velocities in our R — 120, 000 spectrum to a 3a limit of 0.004 A pixel"^. Recently, 
21cm Hi observations have found that high velocity, low column density hydrogen clumps 
have T < 900 K (Hoffman et al. 2004; Richter et al. 2005) and that narrow Call lines are 
detected in most sightlines (Richter et al. 2005). Unfortunately, CaIlAA3935,3970 is not 
covered in either our R = 45, 000 ot R = 120, 000 spectrum of PG 1634 + 706. However, 
at i? = 45, 000, it would not be possible to test if the Call lines arc as narrow as would 
be suggested by the low temperatures from Hi. Thus the study of HVCs is an application 
for which super-high resolution will be particularly valuable. This applies to HVCs around 
galaxies at high redshift as well as to Milky Way HVCs. 



^We adopt the "standard" definition of the local standard of rest (Kerr & Lynden-Bell 1986), in which 
the Sun is moving in the direction a = 18h, S = +30° (epoch 1900) at 20 km s~^. With this convention, the 
conversion to heliocentric velocity is: Vheiio = Vlsr + 14.770 km s~^. 
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4. CONCLUSION 

The ability to routinely obtain quasar spectra with resolutions of > 100, 000 is not 
likely to provide as large a change in the field of quasar absorption lines as the previous 
advance to = 45, 000, made possible with 8-meter class telescopes. However, the ability 
to probe structure on scales of ~ 1 km s~^ will certainly be valuable in select applications and 
in detailed studies of particular absorbers. Examples of such applications include probing 
the physical conditions of cold regions in the ISM of galaxies, particularly in DLAs and 
other strong Mgll absorbers. Also, it appears that high velocity clouds, both around the 
Milky Way and at high rcdshift, may have a cold phase. Finally, surveys extending to small 
equivalent widths will be aided by the enhanced sensitivity of higher resolution observations. 

Our comparison of i? = 120, 000 Subaru/HDS observations to R = 45, 000 Keck/HIRES 
observations of various systems toward PG 1634 + 706 show that at the higher resolution, 
more accurate Doppler parameter measurements can be made, and that additional Voigt 
profile components are resolved. However, for the particular Mgll absorbers along the 
PG 1634 + 706 sightline no extremely narrow individual components were measured. It 
was only for the Galactic Nal line in this line of sight that 6 < 1 km s~^ was measured. 
Because such narrow components are expected in a variety of environments, it would be 
valuable to push the limits of 8-meter class telescopes in order to detect them in a variety 
of types of systems toward the brightest quasars. 

We are grateful to Chris Churchill for giving us permission to use the Keck/HIRES 
spectrum for this work. We also wish to thank the referee, John Black, for several insightful 
comments that have increased the scope of this paper. 
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Fig. 1. — Synthetic absorption profiles of varying Doppler parameter b convolved to spectral 
resolutions of i? = 120, 000 and R = 45, 000. The absorption line is redshifted to 2; = 1. The 
tick marks the center of the absorption profile. 
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Fig. 2. — A sample from our synthetically simulated Mgll A2796 A absorption lines with 
noise added. The two set of lines with Doppler parameters 6(Mg"'") = 2.0 km s~^ and 
6(Mg+) = 1.0 km s-^ have a column density of A^(Mg+) = 10^^ cm-^. The dotted line in 
each panel is the error spectrum. The tick marks represent the center of each absorption 
profile. 
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Fig. 3. — The distribution of measured Doppler parameter b from 10, 000 realizations of an 
Mgll A2796 A line with W^^(2796) = 0.03 A. The solid histor grams are measurements from 
well constrained Voigt Profile models that had a{b) < b. For each case, the measurements 
that were affected by pixelization are separately shown as the unfilled histogram. The true 
value of 6-parameter is 2.0 km s~^ 
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Fig. 5. — The distribution of measured Doppler parameter h from 10, 000 realizations of an 
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Fig. 6. — The distribution of measured A'"(Mg+) column density from 10, 000 realizations of 
the Mgll A2796 A line with Wr{2796) = 0.02 A. The histo gram shown in dotted line are 
measurements that yielded poor constraints (cr(6) > b). The true value of column density is 
7V(Mg+) = 10^2 cni-2. 




Fig. 7. — A sample showing closely blended synthetically simulated Mgll 2796 A hnes with 
noise added corresponding to a S/N ratio. The siib-pancls in the left panel correspond to the 
two absorbing components separated in velocity by 10 km s~^. The right panel are sample 
simulated spectra illustrating the case where the two absorbing components are separated 
in velocity by 7 km s^^. The ticks mark the position of the two components. The S/N ratio 
for the two resolutions are estimated for equal exposure times. 
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Fig. 8. — Synthetic spectra involving isotope and hyperfine lines, without adding Poisson 
noise. The Mgll line is redshifted to z = 1. A column density of A^(Mg"^) = 10^^'^ cm~^ 
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isotopic line shift becomes distinguishable at spectral resolution it! > 500000 and also when 
the Doppler width is smaller. 
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Fig. 10a.~ Comparison of PG 1634 + 706 observations using R = 120,000 Subaru/HDS 
and R = 45,000 Keck/HIRES. The details of the observation are Hsted in Table 1. The 
metal lines associated with the absorption systems discussed in this paper are identified and 
labelled in both spectrum. The features in the HDS spectrum marked with 'x' are unreal 
and were produced by bad regions in the CCD. 




Fig. 10b. — Contd. Figure 10a 




Fig. lOc— Contd. Figure 10a 
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Fig. 11. — The z = 0.9902 strong Mgll system as captured by the Keck/HIRES and Sub- 
aru/HDS spectrum of PG 1634 + 706. The left panel shows the low ionization Mgl feature 
associated with this strong system. The thin solid line superimposed on the spectrum are 
the various Voigt Profile components contributing to the absorption feature. The tick marks 
placed above the features represent the center of each component. The dotted hue is the 
error spectrum. The results from the Voigt Profile fit are listed in Table 3 
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Fig. 12. — The distribution of measured b parameter from 10,000 realization of Mgl 2853 A by 
considering the isotopic hues of Magnesium. The unfilled histogram includes measurements 
that yielded poor constraints((T(6) > b). The vertical line shows 6(^^Mg) = 0.75 km s~^, 
which is the true Doppler parameter for that line. The line was chosen to have a column 
density of A^(Mg) = 10^^'^ cm~^ to match with the model predicted by (Ding et al. 2003). 
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Fig. 13. — The two single-cloud weak Mgll systems at redshifts 0.8181 {left panel) and 0.9056 
{right panel) as detected in the Keck/HIRES and Subaru/HDS spectrum of PG 1634 + 706. 
The dotted line is the error spectrum. The Voigt Profile fit results are hsted in Table 4 
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Fig. 14. — The multiple-cloud weak Mgll system ai z = 1.0414 as detected in the 
Keck/HIRES and Subaru/HDS spectrum of PG 1634+706. The thin solid line superimposed 
on the spectrum are the various Voigt Profile components contributing to the absorption fea- 
ture. The tick marks placed above the features represent the center of each component in 
the absorption profile. The dotted line represents the error spectrum. The Voigt Profile 
fitting results are listed in Table 4 
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Fig. 15. — The z — 1.0414 multiple-cloud weak Mgll system and the associated subsystem 
in the spectrum of PG 1634 + 706. The dotted line represents the error spectrum. The thin 
solid lines mark the respective positions for the AA 2796 A and 2803 A lines that form the 
doublet pair. The two absorption systems are separated by ~ 150 km s^^. In the HIRES 
spectrum, for the bluer subsystem, Mgll A2796 transition was not covered due to an cchellc 
order break. However, the weaker pair of the doublet, Mgll A2803, was not detected down 
to an equivalent-width limit of 3a. However, there is a clear detection of this weak subsytem 
in the HDS spectrum. 
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Fig. 16. — Galactic Nal detections in the HDS spectrum of PG 1634 + 706. The superhigh 
resolution spectrum allowed measurements of 6- value ~ 1 km s^^ Superimposed in thin 
hnes are the various Voigt Profile components contributing to the absorption. The vertical 
ticks mark the center of each component. The thin line in the Nal 5898 panel shows the 
cumulative Voigt Profile by combining the various components. The fit results are listed in 
Table 5. The Mgll absorption from the HVC Complex C along the line of sight to the quasar 
is also shown. To guide the eye, the expected location of the Nal doublet absorption from 
the HVC is marked by the vertical line. 
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Table 1. Results from simulations 



i?/1000 


S/N ratio 




6(Mg+) = 2.0 km 






6(Mg+) = 1.0 km 






pixel" ^ 


6(Mg+)* 


log N{Mg+)* 


% Models t 


KMg+)* 


log 7V(Mg+)* 


% Models t 


120 


35 


2 02+°-^^ 

-'•Ll^_0.14 


19 r^+0.02 
■'■^•'-'-0.02 


99.99 


-1 f.7+0.16 


11.97tro! 


90.21 


45 


80 


9 1 c+0.43 
^•■'-"3-0.48 


11 qq+0.04 
ii.yy_Q Q3 


94.33 


-, 09+0.44 
■'■•°^-0.33 


11.84iro^ 


44.00 


120 


55 


9 no+0.10 
-^•"-"^-0.09 


i2.otro^ 


99.98 




11.97trol 


97.80 


45 


128 


9 -17+0.26 
-0.32 


ii.Qstn? 


98.44 


i-56t°;f7 


11.85tro^ 


51.86 


120 


70 




i2.otroi 


99.98 


i-UO_o.io 


ll-97tro^ 


98.98 


45 


160 


9 -17+0.21 
'^•-'■'-0.25 


ii.Qstn? 


98.03 


-1 cc+0.33 

1.00_o.28 


11.86tro^ 


51.98 



Note. — * Median h and A'' values. 
''' - Out of 10,000 realizations, the percentage of Voigt Profile fit models that gave a reliable measurement 
< h). 
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Table 2. PG 1634 + 706 Observations; ^qso=1-337, my=14.9. 



Instrument 


Resolution 


Av 


Datc(s) of 


T 

J- exp 


Coverage 


PI 




A/AA 


km s~^ 


Observation 


seconds 


A A 




Keck/HIRES 


45, 000 


6.6 


1994 Jul 04 


2700 


3723 - 6185 


Churchill 




45, 000 


6.6 


1994 Jul 05 


5400 


3723 - 6185 


Churchill 


Subaru/HDS 


120, 000 


2.5 


2005 Jun 06 


3600 


4830 - 7560 


Misawa 
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Table 3. VoiGT Profile Fit Results for z = 0.9902 Strong Mgii System 



Instn iment 


V 




h 






km 


atoms cm~^ 


km s~ 


-1 




Mgl 


Fit Results 






HIRES 


-22.6 


11.35 ± 0.04 


6.52 ± 


0.71 




5.2 


11.37 ± 0.06 


13.52 ± 


2.44 




25.0 


10.65 ± 0.18 


3.2 ± 


2.52 


HDS 


-20.1 


11.26 ± 0.03 


5.87 ± 


0.46 




1.1 


11.28 ± 0.04 


10.55 ± 


1.28 




24.5 


10.51 ± 0.10 


5.87 ± 


0.94 




Mgii Fit Results 






HIRES 


-22.1 


13.04 ± 0.01 


7.72 ± 


0.09 




1.1 


12.82 ± 0.01 


4.11 ± 


0.12 




8.3 


13.30 ± 0.01 


22.29 ± 


0.21 


HDS 


-20.2 


13.05 ± 0.01 


7.57± 


0.15 




-8.2 


12.44 ± 0.06 


4.46 ± 


0.48 




3.8 


12.90 ± 0.22 


6.50 ± 


1.00 




16.8 


12.93 ± 0.34 


11.03 ± 


8.00 




30.1 


12.34 ± 0.55 


6.11 ± 


1.91 




45.4 


11.34 ± 0.09 


3.37 ± 


0.69 
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Table 4. VoiGT Profile Fit Results for Weak Mgii Systems 



Tn stn 1 m pnt 


V 


log [iV(Mg+)] 


6(Mg+) 




km s"^ 


atoms cm~^ 


km s~^ 


z 


= 0.8181 


single-cloud system 


HIRES 





12.04 ± 0.03 


2.14 ± 0.40 


HDS 





11.97 ± 0.02 


2.61 ± 0.13 


z 


= 0.9056 single-cloud system 


HIRES 





12.47 ± 0.01 


2.77 ± 0.12 


HDS 





12.42 ± 0.01 


3.14 ± 0.07 




: 1.0414 multiple-cloud syst^ 


em 


HIRES 


-38.2 


11.49 ± 0.08 


6.53 ± 1.62 




-15.2 


11.90 ± 0.06 


10.81 ± 2.03 







12.02 ± 0.03 


3.12 ± 0.38 




15.7 


11.68 ± 0.04 


5.95 ± 0.83 


HDS 


-36.4 


11.17 ± 0.23 


3.98 ± 2.78 




-29.6 


10.98 ± 0.27 


2.23 ± 1.46 




-17.3 


11.38 ± 0.14 


5.19 ± 1.99 




-10.1 


11.05 ± 0.25 


2.35 ± 1.08 




1.9 


12.02 ± 0.01 


3.82 ± 0.13 




18.0 


11.48 ± 0.04 


2.86 ± 0.34 
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Table 5. VoiGT Profile Measurement of Galactic Nai absorption in the 

HDS SPECTRUM OF PG 1634 + 706 



vlsr 


log [iV(Na)] 


b (Na) 


km 


atoms cm~^ 


km s~^ 


-1.1 


10.56 ± 0.10 


1.68 ± 0.82 


-10.0 


11.51 ± 0.02 


1.63 ± 0.16 


-UA 


10.67 ± 0.10 


0.89 ± 0.88 


-28.6 


10.90 ± 0.11 


2.80 ± 1.02 


-33.2 


11.14 ± 0.05 


0.90 ± 0.33 


-49.9 


11.96 ± 0.34 


4.64 ± 0.28 


-51.7 


14.21 ± 1.63 


1.02 ± 0.35 


-59.6 


11.26 ± 0.05 


0.86 ± 0.21 



